Spin density distribution in a partially magnetized organic quantum magnet 



A. Zheludev,i'El V. O. Garlea/ S. Nishihara,^'^ Y. Hosokoshi,^'^ A. Cousson,^ A. Gukasov,"^ and K. Inoue^ 

^Neutron Scattering Sciences Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831-6393, USA. 
^Department of Physical Science, Osaka Prefecture University, ,Osaka 599-8531, Japan. 
^Institute for Nano fabrication Research, Osaka Prefecture University, Osaka 599-8531, Japan. 
^Laboratoire Leon Brillouin, CEA-CNRS Saclay, France. 
^Department of Chemistry, Hiroshima University, Hiroshima 739-8526, Japan. 

(Dated: February 6, 2008) 

Polarized neutron diffraction experiments on an organic magnetic material reveal a highly skewed 
distribution of spin density within the magnetic molecular unit. The very large magnitude of 
the observed effect is due to quantum spin fluctuations. The data are in quantitative agreement 
with direct diagonalization results for a model spin Hamiltonian, and provide insight on the actual 
microscopic origin of the relevant exchange interactions. 



I. INTRODUCTION 

Static spin correlations in low-dimensional and molec- 
ular magnets are often severely affected by zero-point 
quantum spin fluctuations. The simplest and most ex- 
treme example is that of an antiferromagnetic (AF) spin 
dimer, where the spin density distribution S(r) is strictly 
zero in the ground state, even in the presence of a small 
aligning magnetic field. In partially magnetized states 
of more complex systems one can expect to find exotic 
non-trivial spin densities that too are strongly modified 
by quantum fluctuations. Experimental studies of these 
correlations can be revealing of the underlying physics, 
and help determine or validate theoretical models used to 
describe such materials. Below we report a direct mea- 
surement of field-induced magnetization densities in the 
novel organic molecular magnet 2-[2',6'-difluoro-4'-(A^- 
tert-butyl-A^-oxyamino)phenyl]-4,4,5,5 -tetramethyl-4,5- 
dihydro-li^-imidazol-l-oxyl, F2PNNNO for short. 

This compound is a prototypical spin-tetramer 
system^ii^ Its molecular building block (Fig. la) con- 
tains only s- and p- elements, but is nevertheless mag- 
netic, thanks to two unpaired electrons that reside in tt* 
antibonding molecular orbitals of the nitronyl nitroxide 
(NN) and the tert-huiy\ nitroxide (tBuNO) groups, re- 
spectively. The material crystallizes in an orthorhombic 
structure, a = 19.86 A, 6 = 14.01 A and c = 13.48 A, 
space group Fbca- In the crystal, F2PNNNO molecules 
are arranged in pairs, so that their tBuNO groups are 
close enough for partial orbital overlap, enabling inter- 
molecular magnetic interactions. The result is a two- 
molecule unit containing four interacting spins (Fig. lb). 
A model Heisenberg Hamiltonian for these S = 1/2- 
tetramers was proposed based on bulk susceptibility and 
high field magnetization data*^ Intra-molecular exchange 
coupling is ferromagnetic, with Jp ~ 35 meV. Inter- 
molecular interactions are AF in nature, of magnitude 
Jaf ^ 5.8 meV.^ In agreement with experiment, this 
model has a unique non-magnetic ground state with to- 
tal spin ^totai = 0, and a gap in the excitation spectrum. 
The magnetization density in the ground state is strictly 
zero in entire space. 
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FIG. 1: (Color online) (a) Molecular structure of F2PNNNO. 
Hydrogen atoms are not shown. The shaded rectangles rep- 
resent S = 1/2-carrying unpaired electrons distributed over 
the nitronyl nitroxide and tert-huiy\ nitroxide groups, (b) A 
schematic representation of the 4-spin Heisenberg Hamilto- 
nian for a two- molecule F2PNNNO unit. Jf and Jaf are 
ferro- and antiferromagnetic exchange interactions, respec- 
tively. Vertical arrows represent individual spins in the clas- 
sical ground state | TTii)- The actual quantum ground state 
has zero spin density throughout the tetramer. 



A non-trivial spin (magnetization) density is only to 
be found in some of the tetramer's excited states. In the 
presence of an external magnetic field applied along the 
z axis, the one with the lowest energy has a total spin 
5'totai = 1 and a spin projection 6';^, total = +1- We shall 
denote this state as |1,-|-1). By numerically diagonal- 
izing the 4-spin Heisenberg Hamiltonian we find that it 
actually is a linear combination of four "pure" spin wave 
functions: 

|i,+i) = a| TTn) + /3| TnT)-ce| iTTT)-/5| TiTT), (i) 

where a ~ 0.46 and P ^ 0.54. The most striking conse- 
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FIG. 2: (a) Calculated field-dependence of the energy lev- 
els of an F2PNNNO spin-tetramer. Heavy solid lines are the 
lowest energy S = and S — 1 states. Dashed lines are S — 2 
states. More excited states are present above 35 meV. The 
plot is laid over a filled-contour plot of the thermal popula- 
tion function exp(— /ia;/T) for T = 10 K (0.1 contour step), 
(b) Calculated temperature dependence of the sum (dashed 
line) and difference (solid line) of the spin populations of the 
nitronyl nitroxide and tert-h\xiy\ nitroxide groups. The aroow 
indicates the experimental temperature. 



quence of the quantum-mechanical nature of this state is 
a skewed spin density distribution ^'^(r). The local spin 
populations of the NN groups are expected to be equal, 
but different from those of the tBuNO groups. Their ra- 
tio R is given hy R = o? j 0^ ^ 1.39. The central purpose 
of this work is an experimental detection of this effect. 



II. EXPERIMENTAL 

To observe this phenomenon, one must first prepare 
the tetramer in its first excited state. One strategy 
is to substantially increase the external field. Due to 
Zeeman effect, the energy of the excited state will de- 
crease, and eventually reach zero at some critical field 
Re At this point it will become the new ground state, for 
which the spin density distribution can be measured. In 
F2PNNNO, due to residual inter-tetramer interactions, 
the transition at Rc is spread out between i^i = 9 T and 
R2 = 15 T.^ While it is certainly possible to perform 
experiments at > i^i, the equipment available for the 
present study was limited to fields up to 7 T. For this rea- 
son we used a slightly different approach. First, a high 
field was used to lower the energy of the |1, +1) state as 
much as possible. The data were then taken at an ele- 
vated temperature of T = 10 K that made this state par- 
tially populated due to thermal fluctuations. Of course, 
states |1, 0) and |1, — 1), as well as other higher- spin states 
got thermally excited as well. However, at T = 10 K ther- 
mal populations of higher-spin states are negligible. This 
is illustrated in Fig. [2l^a) that shows the field dependence 



of tetramer energy levels calculated using exact numeri- 
cal diagonalization of the Heisenberg Hamiltonian. The 
plot is laid over a shaded contour plot of the thermal 
population function ex.-p{—huj/T). At = 7 T all but 
the lowest S = 1 energy levels are outside the populated 
region. The |1, 0) state plays no role, as it has ^'^(r) = 0. 
The spin density distribution in |1, — 1) is exactly the re- 
verse of that for |1, +1), and does not affect the imbalance 
between the NN and tBuNO groups. The only adverse 
effect of the finite-T approach is a reduction of the total 
magnetization of the tetramer, that ultimately reduces 
signal strength in any spin density measurement. The 
actual value of T = 10 K was selected to optimize both 
the total tetramer magnetization and the predicted pop- 
ulation difference between the NN and tBuNO groups at 
R = 7 T (Fig.[2]3). A full thermodynamic calculation for 
the 4-spin Hamiltonian predicts R = 1.32 for these con- 
ditions, just slightly less than the ideal value R = 1.39 
for a tetramer purely in the |1, +1) state. 

Measuring the distribution of about 0.5 /ie over two 
molecules with 48 atoms each with angstrom resolution 
is a formidable experimental challenge. It can only be 
met by polarized neutron diffraction.^ This technique 
achieves great sensitivity by exploiting the interference 
between magnetic and nuclear scattering of neutrons in 
the crystal. For F2PNNNO the data were taken at the 
5C1 and 6T2 lifting counter diffractometers installed at 
the Orphe reactor at LLB, using 0.841 A and 1.4 A neu- 
trons, respectively. Beam polarizations of 91% or 97% 
were achieved using Heussler-alloy monochromator and 
supermirror bender. A 20 mg F2PNNNO single crystal 
sample was mounted consecutively with the a, b and c 
axes parallel to the field direction. Sample environment 
was a split-coil cryomagnet. Overall, 70 independent flip- 
ping ratios were measured in magnetic fields R = 7 T 
and = 4 T at T=10 K, typically counting 2 hours 
per reflection on the 5C1 and 30 minutes on the 6T2 
diffractometers, respectively. Extracting the correspond- 
ing spatial Fourier components of ^'^(r) from these data 
required knowledge of the low-temperature crystal struc- 
ture. The latter was measured in a single crystal unpo- 
larized neutron diffraction experiment. 3887 independent 
Bragg intensities were measured for a 5 mg single crystal 
sample on the 5C2 4-circle diffractometer at Orphe us- 
ing 0.832 neutrons. Sample environment was a gas flow 
cryostat, and the data were taken at T=50 K. The crys- 
tal structure was refined assuming isotropic vibrational 
parameters for hydrogen atoms, and general anisotropic 
ones for all other. The resulting least-squares R-factor 
was 0.089. 



III. RESULTS 

Inverting the Fourier transform to reconstruct the real- 
space Sz{t) function is far from straightforward, and 
prompted us to apply several complimentary approaches. 
One such tool was the maximum entropy (ME) method.^ 
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TABLE I: Experimental atomic spin populations (fiB units) 
obtained using the AOE reconstruction method. 
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FIG. 3: (Color online) Experimental spin density distribution 
in a F2PNNNO spin-tetramer at T = 10 K, as reconstructed 
using the Maximum Entropy method, (a) and (b) are pro- 
jections onto the (a, b) and (6, c) (crystallographic planes, re- 
spectively. Overlaid are skeletal representations of F2PNNNO 
molecules as they are positioned within the crystallographic 
unit cell. Areas outlined with thick dashed lines were used to 
estimate the spin populations on the nitronyl nitroxide and 
tert-huty\ nitroxide groups (see text). 



The procedure is model-independent: it uses only the ex- 
perimental structure factors and crystal symmetry as in- 
put, and does not rely on any additional information or 
assumptions. It is known to be particularly effective at 
reconstructing 2D projections of Szi^) onto planes close 
to the principal scattering plane of the diffractometer.^ 
Results of such reconstructions for F2PNNNO are shown 
in Fig. 2. Despite the limited experimental spatial reso- 
lution, one immediately sees that the spin density is pri- 
marily localized around the N and O atoms. Already 
at this stage it is possible to get a crude estimate of 
the redistribution effect. Integrating over the empiri- 
cally chose regions outlined in Fig. 2a, we get R = 1.22. 
However, the actual NN/tBuNO spin population imbal- 
ance is likely to be more pronounced than suggested by 
ME. For all its advantages, the algorithm is known to 
systematically bias the answer towards a more uniform 



distribution. 

An alternative reconstruction method known as atomic 
orbital expansion (AOE)^ lacks the benefit of being 
model-independent, but is free of such a bias and is 
better quantified. It involves refining a parameterized 
model for Sz (r) to best-fit the experimental Fourier data. 
Its application to F2PNNNO was founded on previous 
experiments^iii^i^ and first-principle calculations^ for re- 
lated nitroxides. The NN spin density was described in 
terms of five atomic populations. It was assumed to be 
concentrated in pz Slater- type atomic orbitals of the O, 
N and apical C atoms, the z axis chosen perpendicular to 
the corresponding N-O-C planes. Two more parameters 
were used to quantify the spin density delocalized over 
the Pz' orbitals of the tBuNO N and O atoms, with the 
z' axis perpendicular to the tBuNO 0-N-C plane. One 
additional parameter was used for the sign-alternating 
spin density induced in the phenyl ring by virtue of the 
spin polarization effect ^^^^ It was assumed to be contained 
in pz" orbitals of the phenyl's C atoms {z'' is oriented 
perpendicular to the phenyl plane). The final three pa- 
rameters were the radial exponents of Slater-type atomic 
orbitals for the N, O and C atom types. This model 
yields an excellent least-squares fit to the data collected 
at = 7 T and = 4 T, with = 1-09 and = 1-20, 
respectively. Figure 2 is an isosurface representation of 
the resulting 3-dimensional spin density distribution in 
the tetramer at = 7 T. Individual atomic spin popu- 
lations obtained in the refinement are listed in Table [H 

A very good measure of the AOE's reliability is its re- 
sult for the total tetramer magnetization: m — 0.48(2)/iB 
and m = 0.28(2)/iB, for = 7 T and = 4 T, re- 
spectively, at T = 10 K. These values are consistent 
with existing bulk susceptibility data, and agree well with 
a thermodynamic quantum-mechanical calculation for a 
single tetramer: m = 0.59/iB and m = 0.32/iB, respec- 
tively. With this assurance of the validity of our ap- 
proach, we can finally obtain experimental estimates for 
the imbalance between the NN and tBuNO spin popula- 
tions: R = 1.53(3) and R = 1.51(2), for i7 = 7 T and 
i^" = 4 T, based on AOE model refinements. 
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FIG. 4: (Color online) Experimental spin density distribution in a F2PNNNO spin-tetramer at H — 7 T and T = 10 K, 
as reconstructed using orbital model refinement (see text). The isosurfaces are drawn at 1 x lO^^/is/A^ (orange) and — 1 x 
10"Vb/A^ (blue) levels. The axes show cartesian coordinates in Angstroms. 



IV. DISCUSSION 

A value R > 1 signifies a reduction of uniform in- 
duced magnetization around an antiferromagnetic bond 
in the tetramer, and is to be expected. What is im- 
portant though, is that this effect is hugely magnified 
by quantum correlations. In a classical magnet with 
similar exchange constants, all spins would align them- 
selves in the (x, y) plane and tilt slightly in the field {z) 
direction. It is easy to show that the resulting imbal- 
ance in Sz{t) would be an order of magnitude smaller: 

i^classical ^ 1 + | Jaf|/4| Jf| = 1-04. 

The quantitative agreement between experiment and 
direct diagonalization calculations is an important micro- 
scopic validation of the model Heisenberg Hamiltonian 
that was initially hypothesized based on bulk measure- 
ments alone ^2 In fact, our measurements can be viewed 
as a direct experimental determination of Jf/Jaf in 



F2PNNNO. In addition, these experiments help under- 
stand the microscopic interactions within the F2PNNNO 
molecule. In Fig. 3, note the negative density in the 
vicinity of the apical carbon atom of the NN group. This 
large negative spin population^ plays a key role in the 
ferromagnetic intra- molecular coupling Jp. It is a part 
of a sign-alternating spin density wave that propagates 
across the phenyl ring and connects the positively pop- 
ulated N sites of the NN and tBuNO fragments over a 
large distance. This density- wave mechanism is analo- 
gous to Ruderman-Kittel-Kasuya-Yosida interactions in 
metals^iii 



V. CONCLUDING REMARKS 

Field-induced spin distributions in quantum magnets 
are strongly affected by quantum correlations. They can 
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be directly probed by polarized neutron diffraction and 
carry valuable information on the system. A very promis- 
ing avenue for future work are similar experiments on 
F2PNNNO conducted at low temperature in the magne- 
tization plateau phase H > Hc2 and in the inter-plateau 
region Hd < H < Hc2' In these regimes the system is 
expected to be ordered in three dimensions due to inter- 
tetramer interactions. How does transverse long-range 
order influence the distribution of Sz , and is in any differ- 
ent from that in effectively isolated tetramers, as studied 
in this work? 
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